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U ■ ABSTRACT 



We present results from the Surface Brightness Fluctuation (SBF) Survey for the 
distances to 300 early-type galaxies, of which approximately half are ellipticals. A mod- 
est change in the zero point of the SBF relation, derived by using Cepheid distances to 
spirals with SBF measurements, yields a Hubble constant Hq = 77±4±7 kms -1 Mpc -1 , 
somewhat larger than the HST Key Project result. We discuss how this difference arises 
from a different choice of zero point, a larger sample of galaxies, and a different model 
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for large scale flows. The zero point of the SBF relation is the largest source of uncer- 
tainty, and the SBF value for Ho is subject to all the systematic uncertainties of the 
KP zero point, including a 5% decrease if a metallicity correction is adopted. 

To analyze local and large-scale flows — departures from smooth Hubble flow - 
we use a parametric model for the distribution function of mean velocity and velocity 
dispersion at each point in space. These models include a uniform thermal velocity 
dispersion and spherical attractors whose position, amplitude, and radial shape are free 
to vary. Our modeling procedure performs a maximum likelihood fit of the model to 
the observations. 

Our models rule out a uniform Hubble flow as an acceptable fit to the data. Inclusion 
of two attractors, one of which having a best fit location coincident with the Virgo 
cluster and the other having a fit location slightly beyond the Centaurus clusters (which 
we refer to by convention as the Great Attractor) reduces X 2 /N from 2.1 to 1.1. The 
fits to these attractors both have radial profiles such that v ~ r -1 (i.e. isothermal) 
over a range of overdensity between about 10 and 1, but fall off more steeply at larger 
radius. The best fit value for the small scale, cosmic thermal velocity is 180 ± 14 km s -1 . 

The quality of the fit can be further improved by the addition of a quadrupole 
correction to the Hubble flow. The dipole velocity offset from the CMB frame for 
the volume we survey (amplitude ~ 150 km s -1 ) and the quadrupole may be genuine 
(though weak) manifestations of more distant density fluctuations, but we find evidence 
that they are more likely due to the inadequacy of spherical models to describe the 
density profile of the attractors. Our models can account almost perfectly for the CMB 
motion of the Local Group as arising from the attractors within our survey volume 
(R < 3000 km s -1 ); in other words, our sampled volume is, in a mass averaged sense, 
essentially at rest with respect to the CMB. This contradicts claims of large amplitude 
flows in much larger volumes that include our sample. 

Our best-fitting model, which uses attenuated isothermal mass distributions for the 
two attractors, has enclosed mass overdensities at the Local Group of 7 x 10 14 M 
for the Virgo Attractor and 8 x 10 15 M Q for the Great Attractor. Without recourse 
to information about the overdensities of these attractors with respect to the cosmic 
mean we cannot provide a good constraint on Q M , but our data do give us accurate 
measurements of SQ% 3 = 0.33 for the Virgo Attractor and 5Q% 3 = 0.27 for the Great 
Attractor. 

Subject headings: galaxies: distances and redshifts — galaxies: clusters: individual 
(Virgo, Centaurus) - cosmology: distance scale - cosmology: large-scale structure of 
universe 
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1. Introduction 

There is now convincing evidence that the evo- 
lution of large-scale structure is driven by inter- 
galactic dark matter. Understanding the nature 
of this dark matter remains one of the major 
unsolved problems in astronomy. It is for this 
reason that mapping the overall mass distribu- 
tion, independent of the distribution of luminous 
matter in galaxies, must be considered a funda- 
mental endeavor in cosmological research. Unfor- 
tunately, few tools are available for the purpose: 
weak gravitational lensing can trace, statistically, 
the presence of intervening dark matter at large 
distances, and X-ray emitting gas can be used to 
map the gravitational potential of the dark mat- 
ter that binds rich clusters. However, for a point- 
by-point comparison of the density of dark and 
luminous matter, measuring peculiar velocities — 
departures from pure Hubble flow caused by an 
uneven distribution of dark matter — is the only 
effective method available. 

Using galaxies as test particles to sample the 
local gravitational field requires accurate knowl- 
edge of their positions in space, because galaxy 
motions are dominated by the expansion of the 
universe for all but the closest objects. Since 
most methods deliver an accuracy that (at best) 
scales proportionately with the distance, there is 
a premium on more accurate distance measure- 
ments driven by the desire to measure peculiar 
velocities over a large volume of space. The pe- 
culiar velocity that is known best, by far, is the 
370 km s~ x motion of our Sun with respect to the 
cosmic microwave background (CMB), known 
with exquisite accuracy from the COBE mea- 
surement of the CMB dipole anisotropy, which is 
then converted with less certainty to the 600 km s" 
motion of our Local Group of galaxies with re- 
spect to the CMB. A study of peculiar veloci- 
ties for other galaxies in the local universe can 
be thought of as an exercise with two comple- 
mentary aims: (1) to map the dark matter dis- 
tribution in the local universe and compare it 



with the local galaxy distribution (in the pro- 
cess, measuring a representative value of the 
cosmic density parameter £Im) an d (2) to ac- 
count for our Galaxy's motion as the consequence 
of anisotropics of the dark matter distribution. 
Such knowledge has application in interpreting 
large-scale structure inferred from the galaxy dis- 
tribution alone, and to measuring the degree and 
nature of "bias" that may exist between the dark 
and luminous matter distributions, so important 
for testing theoretical models of structure forma- 
tion. In addition, it is likely that a detailed com- 
parison of the dark- matter /baryon/galaxy distri- 
butions on the scale of the correlation length for 
galaxies will shed light on the processes of galaxy 
formation and the role of dark matter in that pro- 
cess. 

The observation of Rubin et al. (1976) first 
raised the possibility of a large-scale deforma- 
tion in the local Hubble flow. Unfortunately, 
the method of distance used in this pioneering 
work was too crude to fight its way through 
the Malmquist-like systematic errors that are 
endemic to the field. The first solid detec- 
tions of large-scale flows and inferred dark struc- 
tures began around 1980, including Tonry (1980), 
Schechter (1980), Yahil et al. (1980), Aaronson 
et al. (1980), Tonry & Davis (1981), and Aaron- 
son et al. (1982), using spiral and elliptical galax- 
ies to map the Local Supercluster. All studies 
detected the infall pattern that was anticipated 
for a sizeable overdensity roughly centered on 
the Virgo cluster, but with an amplitude at the 
position of the Local Group that ranged from 
125 kms" 1 (Yahil ct al. 1980) to 480 kms" 1 
(Aaronson et al. 1980). As we shall see, the flow 
within the Virgo supercluster is much more com- 
plicated than these early models permitted, and 
this range of "infall velocity" is not surprising. 
On thing that was clear, however, was that in 
both amplitude and direction the pull of the Lo- 
cal Supercluster was insufficient to be the sole 
cause of the CMB dipole anisotropy, i. e., the 
motion of the Milky Way or Local Group with 
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respect to the CMB. An additional pull, in the di- 
rection of the Centaurus supercluster, implicated 
as well by the high galaxy density in this direc- 
tion, was suggested by Shaya (1984), Tammann 
& Sandage (1985), and Aaronson et al. (1986). 
Of particular relevance for the discussion of our 
own results is the suggestion by Lilje, Yahil, & 
Jones (1986) that the effect of this more distant 
mass can be seen as a quadrupole term in the 
peculiar velocity field of the Local Supercluster. 

Based on a refined distance estimator for el- 
liptical galaxies, actually a projection of the fun- 
damental plane (Dressier et al. 1987, Djorgovski 
& Davis 1987), Lynden-Bell et al. (1988) used 
a sample of over 400 elliptical galaxies to map 
a much large volume of space. This study con- 
firmed the infall pattern toward Virgo, but cred- 
ited a considerable fraction of the peculiar veloc- 
ity measured for local galaxies to the quadrupole 
of a more distant mass concentration. Surpris- 
ingly, however, the Lynden-Bell et al. study did 
not find the Centaurus cluster to be the center of 
a more distant infall pattern, but put the center 
at a distance approximately 50% greater, with 
Centaurus and its associated groups themselves 
falling into the more distant attractor at a ve- 
locity of order 1000 kms -1 . The pull on the 
Local Group from this "behind Centaurus" at- 
tractor was reckoned to be of order 500 kms -1 , 
perhaps twice as big as the Virgo pull but from 
a velocity distance of around 4000 kms -1 , over 
three times as distant. The implied order-of- 
magnitude greater mass earned it the nickname 
"Great Attractor" (GA). 

However, doubts have remained about whether 
the GA is indeed a well defined overdensity with 
such a large gravitational influence. While Dressier 
& Faber's (1990a, 1990b) distance measurements 
of additional elliptical and spiral galaxies ap- 
peared to show the distinctive S-shaped pattern 
of infall, including "backside infall," the larger 
spiral sample of Mathewson, Ford, k, Buckhorn 
(1992) did not. Instead, these authors argued for 
a continuing high amplitude flow beyond the dis- 



tance identified as the GA center, perhaps the re- 
sult of a more distant gravitational pull, such as 
the "Giant Attractor" suggested by Scaramella 
et al. (1989) associated with the Shapley concen- 
tration of rich clusters, more than three times as 
distant as the GA. 

Further evidence that the GA's role in pro- 
ducing the local flow pattern might have been 
overestimated is discussed by Courteau et al. 
(1993). Adding a more complete sample of spi- 
rals, particularly in the Perseus-Pisces region, led 
to this study's finding of a considerable ampli- 
tude "bulk flow" — that is, a non-converging flow 
whose source, if gravitational in origin, would be 
on a scale considerably larger than 6000 kms -1 , 
a region encompassing the GA. However, a re- 
cent study of spiral galaxy distances by Gio- 
vanelli et al. (1999) contradicts this result by 
claiming that all of the Local Group's motion 
with respect to the CMB can be accounted for 
by sources within the V = 6000 km s -1 sphere. 
Likewise, the Type la supernovae observations of 
Riess, Press, & Kirshner (1995) suggest a more 
distant frame that is at rest with the CMB. 

This lingering uncertainty in the sources of 
the local flow pattern — that is, the contribu- 
tion of the Local Supercluster, the GA, Perseus- 
Pisces, the Shapley concentration, and structures 
beyond — have been made all the more important 
by studies of Lauer & Postman (1994), Willick 
(1999), and Hudson et al. (1999) who find large, 
but not always consistent bulk flows, over even 
larger volumes of the local universe. Indeed, 
these flows have amplitudes large enough, over 
such large volumes, that they appear inconsistent 
with the small scale anisotropy AT/T ~ 10~ 5 
measured for CMB fluctuations. 

The surface brightness fluctuation (SBF) method 
(Tonry & Schneider 1988) of measuring early- 
type galaxy distances is perhaps the most promis- 
ing means for understanding these issues and 
resolving the disagreements. A recent compre- 
hensive review of the SBF method is given by 
Blakeslee, Ajhar, & Tonry (1999). SBF offers 
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an accuracy several times greater than either 
Tully-Fisher or Fundamental Plane distances; as 
a consequence, its susceptibility to Malmquist- 
like biases is reduced by an order-of-magnitude. 
SBF also includes an implicit correction for vari- 
ations in age and metal abundance, in the form of 
color term, which previous methods for measur- 
ing distances to elliptical galaxies did not include. 
The data we present here will have bearing on 
all the issues raised above, but they are limited 
to relatively nearby galaxies (V < 4000 km s -1 ) 
for all but a handful of observations with the 
Hubble Space Telescope (HST). Nevertheless, we 
hope to demonstrate in this paper the power and 
promise of the method for eventually produc- 
ing a highly accurate map of the local distribu- 
tion of dark matter, out to velocity distances of 
10,000 kms -1 or more, thereby answering the 
question "over what scale does the motion of the 
Local Group with respect to the CMB arise?" 

In this paper we will discuss the characteristics 
of our SBF data set and use the data to construct 
parametric models of the local flow pattern. We 
will show not only that the convergent flows into 
the Local Supercluster and Great Attractor dom- 
inate the departures from smooth Hubble flow 
but also that the data are sufficiently accurate 
to provide other constraints on bulk flow, local 
voids, and other possible attractors. The com- 
bination of these two attractors and a moderate 
thermal component accounts for 90% of the vari- 
ance in our sample. Specifically, evidence for the 
GA comes not only from the large peculiar ve- 
locities in the Centaurus direction but also from 
the clear quadrupole signature in the Local Su- 
percluster flow, whose amplitude and direction 
are consistent with an order-of-magnitude larger 
mass concentration in the Centaurus direction. 
The amplitude and scale of these two flows are 
consistent with Iow-^m cosmologies and with the 
measured small-scale anisotropy of the CMB. 

We extensively explore the covariances be- 
tween various model parameters, such as the at- 
tractor distances and amplitudes, the Hubble pa- 



rameter, and the overall bulk flow. The un- 
certainties will greatly diminish when additional 
HST SBF measurements spanning the GA region 
are available. A future paper will discuss more 
detailed models for the GA and compare our 
results to previous works, including a point-by- 
point comparison of SBF and Lynden-Bell et al. 
distances. In addition, while we have chosen on 
this first examination of our data to use only 
the peculiar velocity information inherent in the 
SBF distances and to completely ignore the dis- 
tribution of galaxies, we intend in future papers 
to compare the SBF peculiar velocities in detail 
with the galaxy density field. 

2. The Data 

2.1. The SBF Distance Survey 

The first paper of this series (Tonry et al. 
1997, hereafter SBF-I) described the galaxy sam- 
ple and observations of the SBF survey in detail 
and gave a new calibration of the method. That 
calibration was derived by comparing averaged 
SBF measurements within 7 galaxy groups to 
Cepheid distances to spirals purportedly within 
the same groups. It gave an absolute magnitude 
= —1.74 at our fiducial early- type galaxy 
color of (V—I)q = 1.15. The alternative, di- 
rect calibration from SBF measurements for 5 
Cepheid-bearing galaxies gave a brighter zero 
point of M°j = —1.82. As these two zero points 
agreed within the uncertainty, and because of the 
potential for systematic differences arising from 
the relative difficulty of SBF measurements in 
spiral bulges, we felt more secure in adopting the 
group calibration. 

For a number of reasons it has become nec- 
essary to revise the calibration yet again. First 
of all, we wish to incorporate all the latest HST 
Cepheid distances, and we adopt the values as 
tabulated by Ferrarese et al. (1999). These au- 
thors have made a number of small revisions 
to the earlier results for purposes of homogene- 
ity. A bigger change is that we have switched 
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to the new extinctions of Schlegel, Finkbeiner, 
& Davis (1998, hereafter SFD), which are deter- 
mined from the COBE/DIRBE and IRAS/ISSA 
dust maps. These have better spatial resolution, 
uniformity, and accuracy than the old estimates 
derived from Hi maps (Burstein & Heiles 1984, 
hereafter BH). While there is a good correla- 
tion between the two sets of reddening estimates, 
there is also a significant zero-point offset. SFD 
find a mean extinction that is greater by 0.02 mag 
in E(B—V) than that of Burstein and Heiles. 
Therefore when we correct our SBF magnitudes 
and colors according to these new extinctions, 
our Mi zero point should change. 

Increasing the extinction estimate to any given 
galaxy (while keeping a fixed calibration) causes 
the SBF distance to that galaxy to increase. This 
is due to the sensitivity of Mj to (V—I) color; 
if the extinction is changed by 5E(B—V), the 
distance modulus changes by 

5{m-M) = [4.5 (A v -Ai) - Ai) x SE{B-V) , 

r (1) 

where Ay and Ai are the ratios of the total V and 
/ band extinctions to the selective E(B—V) ex- 
tinction. For the average increase of 8E(B— V) ~ 
+0.02 in going from BH to SFD extinctions, we 
find a mean 5{m— M) ~ +0.1 mag. However, if 
all the calibrator galaxies underwent this same 
average shift in extinction, then the zero point 
we would derive through comparison with the 
Cepheids would change by 6Mi ~ +0.1 mag, and 
our overall distance scale would not change. This 
is basically the case for the group-derived calibra- 
tion, which now gives M® = — 1.61 ± 0.03 from a 
total of 6 groups. However, the direct Cepheid- 
galaxy calibration did not change, since the mean 
BH and SFD extinctions towards these spirals 
are the same. For M31, which is excluded from 
the SFD map, we follow Bianchi et al. (1996) 
and Ferrarese et al. (1999) in adopting the BH 
value of E(B—V) = 0.08, and then we use the 
SFD extinction ratios for consistency. Thus, the 
galaxy calibration still gives a weighted average 
of Mi = —1.80 + 0.08, from 6 individual compar- 



isons now, and a median of Mj = —1.74. The 
group and galaxy calibrations have diverged from 
a 0.08 mag difference in SBF-I to a 0.13-0.19 mag 
difference with the SFD extinctions. 

In contrast to the approach of SBF-I, we have 
decided to adopt the direct galaxy calibration 
rather than the group calibration. This fairly 
major shift is inspired by new worries that the 
HST Cepheid-bearing spirals may be system- 
atically in the foreground for Virgo and other 
groups (see Ferrarese et al. 1999). Since we are 
not confident about the formal error estimates 
for either the spiral SBF or Cepheid distances, we 
adopt the median rather than the weighted av- 
erage of the spiral zero points. This new calibra- 
tion is within 0.07 mag of the theoretical calibra- 
tion from stellar population models of Worthey 
(1994) (see discussions by SBF-I and Blakeslee 
et al. 1999). The new SBF calibration derived 
with all the latest measurements and the SFD 
extinctions and used in this paper is 

Mi = -1.74 + 4.5 [(V-J)o- 1.15]. (2) 

It is slightly fainter than the SBF calibration of 
Ferrarese et al. for two reasons: they used the ex- 
tinction ratios from Cardelli, Clayton, & Mathis 
(1989) whereas we have used the ones recom- 
mended by SFD, and we have used a median 
rather than an average. Despite the numerical 
coincidence of our new calibration with SBF-I, 
the change due to the increase in average extinc- 
tion [Eq. (||)] means that our distance moduli 
increase on average by ~0.1 mag. We provide 
complete details about our choice of zero point 
in Appendix B. 

Finally, the survey data reductions were not 
yet complete at the time of SBF-I. In particu- 
lar, we have added a large body of new measure- 
ments from the 2.4 m Las Campanas telescope 
that were generally taken under superior condi- 
tions and have enabled us to improve some dis- 
tances. From detailed comparisons of these new 
data to older, poorer quality data, we discovered 
that bad data sometimes yield spurious SBF sig- 
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nals. We have found two simple quantities useful 
in comparing data qualities: pd = psf x vcmb, 
where PSF is the full-width at half-maximum of 
the point spread function in arcseconds and vqmb 
is the CMB velocity in units of 1000 kms" 1 , and 
q = log 2 (#e/m/PD 2 ), where #e/m is the num- 
ber of electrons that would be detected in the 
image from an object of magnitude m (estimated 
from the CMB velocity), and pd 2 is proportional 
to the metric area within a resolution element, 
hence the number of stars for fixed surface bright- 
ness. When pd > 2.7 or q < 0, the observation 
is marked as unreliable and ignored (see § |2.3| ). 
For multiple observations of a given galaxy, we 
adopt a weighted average if the difference in pd 
is less than 0.3 or if each individual observation 
has pd < 1.3; otherwise, we reject the observa- 
tion with the larger pd. 

The entire SBF Survey data set will be pub- 
lished in a forthcoming supplementary paper 
(Tonry et al. 1999, in preparation). 

2.2. HST Data 

We incorporate eight distances from HST SBF 
observations: four from Lauer et al. (1998), 
NGC 4373 from Pahre et al. (1999), and three 
from our Cycle 6 program (Dressier et al. 1999). 
We have adjusted these distances slightly from 
the published ones to agree with the SFD extinc- 
tions and have revised the zero point from Ajhar 
et al. (1997) to agree with our new /-band zero 
point. Although the results presented here do 
not change significantly when the HST data are 
included, they do become more robust to covari- 
ances between, for example, Hq and infall ampli- 
tudes. 

2.3. Biases and Uncertainties 

When we compare these HST distances with 
the difficult ground-based efforts beyond 3000 kms 
or so, we find that there may be a bias in the 
ground-based data that correlates with pd. Typ- 
ical ground-based data for more nearby galax- 
ies agree well with HST measurements (Ajhar 



et al. 1997; see also the discussion by Blakeslee 
et al. 1999). However, at pd = 3 we find that 
dg-b/d-HST ranges between 1 and 0.8, but we do 
not have enough HST observations to be confi- 
dent that this is the onset of a bias in the ground- 
based data. If it truly is a bias, we believe the 
source could have origins including unresolved 
dust and structure, unresolved color gradients, 
and instrumental effects, and we cannot predict 
reliably when it may be present. We therefore 
restrict ourselves to pd < 2.7. When we fit 
our models to subsets of the data with differ- 
ent pd cuts, we find no change in the parameters 
to within the errors for 2.0 < PD < 3.5, which 
gives us confidence that we are not affected by a 
distance-dependent bias. 

As in SBF-I, we try to estimate the accuracy 
of our error estimates by looking at the scatter 
within groups. Since we do not know the ra- 
dial extent of groups, we ask how much cosmic 
scatter in M is necessary if the groups have no 
radial depth at all, and find the answer is 0.10 
mag. Conversely, if we assume that groups have 
the same depth as breadth then the SBF cosmic 
scatter is 0.06 mag. Including this 0.06 mag of 
cosmic scatter in Mi, the quartiles of the distri- 
bution in error in SBF distance modulus are 0.17, 
0.21, and 0.29 mag, so this is a fairly sensitive 
test that the SBF error estimates are accurate. 
(Note that these errors are dominated by the ob- 
servational compromises inherent in doing such a 
large survey on relatively small telescopes; in al- 
most all cases the distances could be significantly 
improved.) 

Because the differential volume increases with 
distance, there will be more galaxies scattered 
down in distance by the errors than are scat- 
tered up. The level of this Malmquist-like bias 
in our sample is difficult to quantify. We have 
the usual problem of not knowing what the true 
distribution of galaxies is, and we may also be 
subject to possible biases and selection effects 
which depend on distance. We try to minimize 
these effects by tailoring our exposure times and 
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seeing conditions to the distance of the galaxy 
we are observing. In addition, we are not actu- 
ally sampling a population that increases rapidly 
with distance (mostly we tend to be seeing early- 
type galaxy in groups). Nevertheless, if we calcu- 
late the e 3 - 5 ( Sd /<f) correction (Lynden-Bell et al. 
1988) for the median error, we get a bias factor 
of 1.038. Including this correction in the models 
described below does not change the results, and 
in fact it slightly increases the \ 2 values. 

The level of bias due to the irreducible cosmic 
scatter in Mj is very small. We preferentially 
measure galaxies to be too close because of their 
cosmic scatter at the 0.3% level (assuming 0.06 
mag of cosmic scatter). Moreover, unlike the 
case for other methods, such as the "forward" 
Tully-Fisher relation, Z? n -cr/fundamental plane, 
or SNIa, that use some of the same observables 
(magnitude, effective radius, surface brightness) 
for the distance measurement as were used in the 
selection, the cosmic scatter in the Mj-(V—I) re- 
lation is not known to correlate with any observ- 
ables measuring galaxy size or luminosity. The 
bottom line is that we can safely ignore bias cor- 
rection for cosmic scatter, and because the dubi- 
ously applicable Malmquist bias factor of ~ 1.04 
yields no measurable improvement, we ignore it. 

The galaxy velocities come from ZCAT (Huchra 
et al. 1992, version dated 1998 November 23), 
and are converted to the CMB frame according to 
Lineweaver et al. (1996), and to the Local Group 
frame as defined by Yahil et al. (1977). (A recent 
study by Courteau & van den Bergh 1999 finds 
almost exactly the same Local Group reference 
frame.) 

There are 336 galaxies in our sample, which is 
reduced to 295 after applying the pd < 2.7 cut. 
We restrict another 10 of these from contribut- 
ing to the model fits, three because we think 
the SBF distance may be in error despite pass- 
ing the pd cut, and seven because their velocities 
are quite discordant with our model velocity dis- 
tribution function. The former cases include a 
messy galaxy at fairly low galactic latitude with 



marginal seeing (NGC 6305), and galaxies with 
extremely blue color where we do not believe our 
calibration is reliable (NGC 5253 and IC 4182). 
The latter cases include three galaxies in Cen-45 
(NGC 4709, NGC 4616, and D 45), NGC 1400 
(the high velocity companion of NGC 1407 in Eri- 
danus), NGC 4150 (a Coma I galaxy with a very 
low velocity), and NGC 4578 and NGC 4419, 
which we think have just passed through Virgo at 
high speed and therefore have extremely unusual 
velocities for their location. 

3. Models of Large Scale Flows 
3.1. Overview 

In this paper we shall limit ourselves to a sim- 
ple, parametric model for the velocity field of 
galaxies. In addition, we try to avoid all group in- 
formation. The assignment of galaxies to groups 
has had a checkered history, and we prefer to deal 
with the virial velocities directly, rather than try- 
ing to average them away. 

There are two methods commonly used for 
constructing a merit function for a peculiar ve- 
locity model V(r). The first is to assemble 

where e 2 = Sr 2 + 5v 2 + a 2 , and 5r is the uncer- 
tainty in position expressed in units of kins" 1 , 
5v is the (small) uncertainty in the measure- 
ment of velocity, and a v is an allowance for ther- 
mal or virial velocity dispersion (cf. Hudson 
et al. 1997). This tacitly assumes that the slope 
dV(r)/dr = H = 1 everywhere, otherwise it is 
not correct to use distance error 5r in the denom- 
inator. This can be very problematic for a model 
which has a significant disturbance to the Hub- 
ble flow, for example in the local supercluster. It 
is possible to "correct" 5r by a spatially varying 
H(r), but then this % 2 is no longer maximum 
likelihood (unless it includes the often neglected 
— 21ncj term) and will tend to bias parameters 
in such a way as to maximize H5r. The other 
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weakness of this approach is that it involves some 
arbitrariness of when objects should be grouped 
and what virial velocity cr v should be allocated 
for them. 

The second common method (e.g., Davis, Strauss, 
& Yahil 1991; Strauss et al. 1992; Shaya, Tully, & 
Pierce 1992) is to invert the velocity model and 
use the observed velocity to provide a model dis- 
tance which can be compared with the observed 
distance. This has the advantage that it com- 
pares the distances and model directly, but it 
has the disadvantage that it cannot cleanly deal 
with virial velocities, necessitating grouping to 
the point that virial velocity is reduced to a neg- 
ligible magnitude compared to the flow velocity 
and distance error. It also has the disadvantage 
that the inversion of V{r) can be multiple valued, 
and is necessarily multiple valued in the regions 
of space where a supercluster flow is most promi- 
nently seen. One is then forced to an arbitrary 
or probabilistic choice of which model distance 
to use. 

What we have chosen to do is to accept the 
fact that there is a distribution function of galaxy 
velocity which varies from place to place, both in 
mean velocity (large scale flow) and in disper- 
sion. Our models consist of a velocity distribu- 
tion function at each location, P(v\r), which we 
take to be Gaussian of mean vo(r) and disper- 
sion cr v (r). A given distance measurement, itself 
a probability distribution, is multiplied by this 
model distribution function and integrated over 
distance, giving a velocity probability distribu- 
tion. The merit function is then the likelihood, 
the product of these velocity distributions evalu- 
ated at the observed velocities. 

The components of our model for the mean 
velocity vo(r) include a Hubble flow of ampli- 
tude Hq, an constant dipole velocity w, possi- 
bly a quadrupole Q with zero trace which acts 
like an anisotropic Hubble constant, a density 
parameter for the universe, and attractors 
which are assumed to be spherical power laws in 
density with a core and cutoff radius. The posi- 



tions of these attractors can be free parameters 
as well as the power-law exponent, the core and 
cutoff radii, and the overall normalization. Our 
model for the velocity dispersion function of 
position, <7 v (r), consists of the quadrature sum 
of an overall thermal velocity, a virial compo- 
nent at the center of each attr actor, and possi- 
bly other (non-attracting) virial components at 
the locations of groups such as the Fornax clus- 
ter. These components have a velocity disper- 
sion which varies spatially as exp(— r 2 /2r^ irial ). 
In principle we have enough distance accuracy 
to constrain both the background thermal com- 
ponent and these virial components (except for 
GA/Centaurus at the limit of our survey) with 
a smoothing length of perhaps 4 Mpc. In prac- 
tice we fit only for the thermal background and 
Virgo velocity dispersions and put in by hand the 
Fornax (235 km s -1 ) and Centaurus (500 km s -1 ) 
dispersions with r v j r j a i = 2 Mpc. Unless other- 
wise stated, we use a thermal velocity dispersion 
of 187kms _1 , which is the best-fit value in our 
most refined models. 

Figure [l] illustrates how this works. The 
model P(v\r) is shown along a line of sight near 
the Virgo cluster. The mean velocity has a Hub- 
ble flow term modified by a peculiar velocity 
which grows approximately as r~ 1 (for an r~ 2 
density distribution), but then rolls off to zero 
at the center of the cluster because of the core 
radius in the density distribution. The velocity 
dispersion at each location consists of the quadra- 
ture sum of the thermal background component 
and a cluster velocity dispersion with an ampli- 
tude of about 650 km s -1 , which declines away 
from the Virgo cluster Gaussian of width 
2 Mpc. Again, these parameters are all variable 
components of the velocity model. 

With three shape parameters (power-law slope 
and core and cutoff radii), our attractor mod- 
els have much radial flexibility. The models can 
emulate anything from a centrally concentrated 
mass to an extended distribution with divergent 
total mass. 
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3000 




Distance (Mpc) 

Fig. 1. — The model velocity distribution func- 
tion P(v\r) is shown for a line of sight which 
passes through the Virgo cluster. The different 
grayscale levels show the 2, 1, 0.5, and 0.2 sigma 
points on the velocity distribution at a given dis- 
tance. A distance observation is shown on the ab- 
scissa as a distance probability function, and the 
resulting velocity probability function is shown 
on the ordinate. This is evaluated at the ob- 
served velocity and forms a term in the likelihood 
product. 

The question of spherical symmetry is harder 
to justify. We note, however, that the galaxy dis- 
tribution in the Virgo supercluster is not com- 
pletely flat, that the potential is much rounder 
than a mass distribution, and that the large 
scale flows we see today arise from the time in- 
tegral of the potential so that a spherical model 
might represent the velocity data quite well even 
though the potential is today significantly flat- 
tened. Moreover, we do allow for the possibil- 
ity of a quadrupole correction to the flow model, 
and this can start to modify a spherical flow into 
a more flattened one. As we shall see, there are 
indications that the Virgo Attractor is not spher- 
ical. 



We convert a mass distribution to a peculiar 
velocity distribution using the usual non-linear 
"nested Friedmann universe" model, where each 
mass shell evolves from the Big Bang according 
to its interior mass density. We do not do this 
computation exactly, however, but rather use the 
Yahil (1985) "p 1 / 4 law" approximation 

Onfall = \h r n% 6 5 {l + 5Y l /\ (4) 

where r is the radius of the shell and 5 is the 
mean mass density interior to the shell in units of 
the background density. Giavalisco et al. (1993) 
examined this approximation with N-body simu- 
lations and found that it is remarkably accurate 
for —0.5 < S < 20, which is the range of interest 
here. We do not expect our models to be valid all 
the way into the virialized core of the attractor 
(5 > 200), so our velocity distribution function 
rolls over to zero in the mean at the centers of 
the attractors with a compensating rise in virial 
velocity dispersion. 

It is common to recast this parametrization in 
terms of a local overdensity and peculiar velocity, 
which then permits £Im to be derived from this 
formula. However, because of our desire to avoid 
use of the galaxy distribution, we lack any a pri- 
ori information about 8, and thus have very lit- 
tle constraint on Qm- Only the non-linear part of 
Eq- @ provides any constraint at all, and clearly 
that is not going to be very reliable. Instead, we 
choose a single value for Q<m and derive the over- 
densities for the model attractors from the fitted 
infall velocities. 

3.2. Details of Models 

The input to our models includes for each 
galaxy the position on the sky, distance and er- 
ror, and velocity. The fitting program reads 
these data and also accepts a variety of model 
components, including cosmology (Hq, Qmi and 
dipole velocity w), quadrupole (five components 
plus origin and cutoff), extended attractors, ther- 
mal components, and compact attractors. The 
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program computes a model velocity distribution 
function for each observed data point, evaluates 
the probability of the observed velocity given 
the velocity probability distribution, forms the 
likelihood sum from the probabilities P: £ = 
In P + constant, and searches parameter space 
for a maximum in C. 

The model velocity field consists of a mean 
velocity at each point and associated Gaussian 
velocity dispersion. The Hubble constant and 
dipole velocity give a contribution to the model 
velocity ^modei at r of w + Hqt. The quadrupole 
contributes a velocity at (x, y, z) of 
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(5) 

The quadrupole matrix is forced to zero trace 
by insisting that Q yy = —Q xx — Qzz, which en- 
sures that the monopole Hq parameter carries 
the net Hubble expansion proportional to dis- 
tance. Unless otherwise specified, the origin of 
the quadrupole, vq, is taken to be the Local 
Group, and the cutoff radius, r qua( j is taken to 
be infinity. 

Extended attractors are modeled by starting 
with spherical density distributions of the form 



p(r) 



(1 + r 3 /r c 3 )7/3 ' 



(6) 



where 7 is the power-law exponent and r c is a 
core radius. This integrates nicely to give an 
enclosed mass function to which we append an 
exponential cutoff and divide by volume to get a 
mean enclosed over density: 



p(r) 



Poe 



—r/r c 



1-7/3 



(1 + 3) W3 



(7) 

Since the central density is finite, the flow veloc- 
ity rolls over at approximately r ~ r c and ap- 
proaches zero as r — ► 0. We ordinarily normalize 
these attractors in terms of a velocity amplitude 
at our location by using the inverse of Yahil's 



p 1 / 4 law to convert such a velocity into an expres- 
sion for pq. The cutoff radius r cut ensures that 
these models have zero net mass, hence do not 
bias Hq high by making a net local modification 
to £Im- For typical values of 7, r cut is approxi- 
mately where the density becomes negative and 
~p(r) starts to decrease back to zero. 

The covariance between the parameters 7, r c , 
and r cut makes it difficult to interpret values for 
any one parameter in isolation. We will strive to 
indicate what is well constrained by the models 
(e.g., mass of attractor, or run of peculiar veloc- 
ity with radius) and what is not. 

As mentioned above, the model velocity dis- 
persion consists of the quadrature sums of back- 
ground and extended thermal components. Each 
component has a velocity dispersion with a cen- 
tral value that falls off spatially according to a 
Gaussian of width r v j r i a i . These thermal compo- 
nents can either be specified explicitly (the cos- 
mic dispersion is taken to be centered on the ori- 
gin with an infinite Gaussian radius; Fornax is 
centered on the location of the Fornax cluster 
and has a core radius of 2 Mpc), or as part of an 
extended attractor. The thermal component of 
an extended attractor is centered on the attrac- 
tor, has a parameter describing its central virial 
velocity dispersion, and shares the core param- 
eter r v i r ; a i with the core radius r c of the mean 
infall. (We expect r viria i to be close to r c , and in 
our models we set the two radii to be the same.) 

It is also possible to specify "Compact Attrac- 
tors" , density distributions which are Gaussians. 
These carry five parameters, their location, their 
Gaussian radius, and their amplitude. The am- 
plitude is given in terms of the total mass and a 
fiducial radius, Rm = 50 Mpc: 



GM tot 1 



J amp 



Rn A Hn 



(8) 
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For a Gaussian with p = po exp(— r 2 /2r 2 ), the 
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enclosed mass is 



M(r) = 4vr/> 



2 VV2ry 



r rexp 




and Mtot = ^Poy f T ■ The expression for M(r) 
is converted to a 5 and then used in Eq. (Q). 
The velocity amplitude can be made negative, 
which corresponds to a negative density void. Of 
course, the void cannot become deeper than the 
mean cosmological density. These compact at- 
tractors do not carry a thermal component. 

We make no attempt to create an overall self- 
consistent density and velocity distribution — these 
components do not interact with one another. 
Having specified the various components, the 
mean velocity is just the sum of the contributions 
of all the components, and the velocity dispersion 
is the quadrature sum of the various contributors. 

Once we have a model, we then go about eval- 
uating how well it matches our observations. For 
Gaussian statistics, we would have that the prob- 
ability for making an observation Vi at a loca- 
tion ri where our model predicts a mean velocity 
V{ji) and Gaussian dispersion a v (ri) is 



1 



2ira v (ri) 



Pi = — — exp -- -— Av . 

2 V cr v {ri) J 

(10) 

The choice of Av is arbitrary but represents how 
closely the model adheres to observation. We 
define the likelihood C = In FJ Pj 



2^V ovfa) 
+ ^-ln\/27rer v (ri) + In Aw. 



(11) 



The first term in — 2C is just the usual definition 
of x 2 , and we can evaluate the goodness of fit 
according to the x 2 P er degree of freedom. If the 
dispersions cr v (rj) do not depend on the param- 
eters, then the second two terms can be ignored, 
since they are constant for a given set of observa- 
tions and do not affect the choice of parameters 
which maximizes the joint probability. 



Our situation is not quite so simple, because 
we believe we have errors which are normally dis- 
tributed in log distance, and although our model 
distribution function has a Gaussian distribution 
of velocity at any point, the mean velocity and 
velocity dispersion vary as a function of posi- 
tion, and hence are not a constant over the range 
where a given galaxy might lie. We deal with 
this for a galaxy of observed modulus fi ± d\x by 
forming a distance probability distribution that 
consists of 11 points between [i — 2dfi to jx + 2d[i 
weighted according to Gaussian statistics. The 
model is evaluated at each of these 11 points, 
providing a Gaussian velocity distribution with 
some mean and dispersion at each point. We 
sum these 11 Gaussians in velocity, weighted by 
their distance probabilities, to form a net prob- 
ability distribution for the velocity we expect to 
see for this galaxy. 

Our likelihood function is then simply formed 
by evaluating this probability at the observed ve- 
locity and summing a negative, normalized like- 
lihood 



-2 J] [in Pi( Vi ) + ln(300 >/2tt)] , (12) 



where the constant term is introduced for Av in 
order to shift the zero point of AT into approxi- 
mate agreement with x 2 ■ Deviations of M about 
its minimum, J\f — N m in = 2(C max — C) are de- 
scribed approximately by a x 2 distribution, so 
we use M to evaluate goodness of fit and confi- 
dence intervals on parameters. We also calculate 
a variance for the velocity distribution function 
we derive for each galaxy and use them to form a 
traditional x 2 ■ This x 2 ls useful for evaluating a 
goodness of fit, but does not correspond to maxi- 
mum likelihood and is not minimized for the best 
fitting model. 

4. Fits of Model to Data 

The model described in the previous section 
carries a lot of parameters, and there can be sig- 
nificant covariance between them, for example 7 
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and r cut , or w x and uga- We will first present 
a sequence of models which motivate our choices 
of mass components, starting with a Hubble flow, 
and ending with two attractors contributing pe- 
culiar velocities. 

We will subsequently look more carefully at 
these covariances and try to show what is well 
constrained by our data and what is not. For ex- 
ample, we cannot independently choose a unique 
7 and r cut , but all our models give nearly the 
same run of peculiar velocity between 6 and 20 
Mpc from the center of the Virgo cluster. Dis- 
tances will often be given in terms of their SGX, 
SGY, and SGZ components in the supergalactic 
coordinate system, defined in the RC3 (de Vau- 
couleurs et al. 1991). 

It is extremely important not to force the 
models to conform to any particular velocity ref- 
erence frame. For example, insisting on the CMB 
frame and using an r~ 2 density model for the 
Great Attractor leads to a very large infall ve- 
locity because of its covariance with w x . We use 
the CMB frame for velocity data, but we always 
fit for an arbitrary dipole vector w as part of our 
models. 

4.1. A First Look 

Figure || shows all of the galaxies for which we 
have SBF distances plotted in a Hubble diagram 
— redshift in the CMB frame as a function of dis- 
tance. The expansion of the universe is apparent, 
and a naive linear fit, without regard for error 
bars or the fact that the abscissa carries greater 
errors, yields a Hubble ratio of 73 kms -1 Mpc -1 . 
Figure || focusses in on the more distant galaxies 
with smaller errors and plots Hubble ratio as a 
function of distance. Local peculiar velocities 
cause scatter in the Hubble ratio at distances of 
40 Mpc, but it seems to settle down fairly well to 
an asymptotic value of 73kms _1 Mpc -1 . How- 
ever, as is seen in Figure ^, there is a lot of scatter 
in the Hubble plot, more than can be explained 
by distance error. 
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Fig. 2. — Recession velocities in the CMB frame 
are plotted as a function of distance. The line 
has a slope of 73kms -1 Mpc -1 ; the very deviant 
high points near 40 Mpc are Cen-45 galaxies. 

The near field Hubble flow is shown in Figure 
U in two different directions. Velocities in the 
Local Group frame are used to avoid the constant 
velocity offset incurred by using the CMB frame 
locally. It is possible to discern the important 
features of the local large-scale flow. The Virgo 
cluster in the lower pane at SGY = +20 Mpc lies 
several hundred kms -1 below the Hubble flow — 
this is often referred to as our "infall velocity" 
towards Virgo. However, this is not the entire 
story, since the Fornax cluster at —20 Mpc also 
has a smaller velocity than the Hubble flow, but 
the quadrupole from a pure Virgo infall should 
cause it to have a larger velocity (although it is 
quite near the quadrupole null at 45°). Likewise, 
the quadrupole from a pure flow towards Virgo 
should cause the Hubble ratio to be smaller in a 
direction perpendicular to Virgo, but in fact we 
find a ratio in the upper panel which is slightly 
larger than the nominal Hubble flow. What we 
are seeing in these local galaxies are the effects 
of two flows, one towards Virgo and the other 
towards Centaurus in the —SGX direction. 
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4.2. H only 

We start by fitting the distances and veloci- 
ties using a model which includes Hq and a pe- 
culiar velocity w, for a total of 4 free parame- 
ters. The sample consists of 285 galaxies. Al- 
lowing 187 kms -1 of thermal velocity but no 
additional virial velocities, we find Hq = 70 
kms -1 Mpc -1 and an overall dipole peculiar ve- 
locity of w = (-330, +180, -80) kms" 1 in the 
CMB frame. The J\f for this model is 513, and 
X 2 per degree of freedom (DOF) for this model 
is x% = 2-08. 

An obvious shortcoming of this model is that 
there is no allowance for virial velocities other 
than 187 kms -1 of background thermal veloc- 
ity. If we add in additional thermal components 
for Fornax and Centaurus (235 and 500 kms -1 , 
which will be used throughout) and Virgo (650 
kms -1 , which is about what we find when we fit 
for it), J\f drops considerably to M = 387 and 
X% = 1.41. In the CMB frame we again get 
Hq = 71 with a negligibly different peculiar ve- 
locity of w = (-330, +200, -90). 

Figure |5] illustrates residuals of velocities af- 
ter removal of this model in the CMB frame. 
In all these vector field plots, the black points 



Fig. 4. — Recession velocities (in Local Group 
frame) are plotted as a function of distance in 
two directions. The upper panel shows galax- 
ies which lie approximately in a 30° (half angle) 
cone in the +SGX direction, 90° away from the 
direction of Virgo. The lower panel shows galax- 
ies which lie within a 15° cone towards Virgo 
(+SGY) and a 45° cone away from Virgo (which 
includes the Fornax cluster). A Hubble ratio of 
73 kms -1 Mpc -1 is also drawn. 

(falling stones) are blueshifted residuals which 
are greater than 1-a, the white points (rising bal- 
loons) are redshifted residuals greater than 1-a, 
and the gray points are residuals which are less 
than 1-a in absolute value. The arrows indicate 
the magnitude of the residual according to the 
Hubble flow, i.e. an arrow of length 5Mpc rep- 
resents a velocity residual of 5Mpc x ffokms -1 . 
All the points are projected onto the supergalac- 
tic x— y plane, and the region where the galactic 
plane cuts through is darker gray. When galaxies 
are members of a group we plot only the group 
residual, so as to prevent the plot from becom- 
ing too overcrowded. The main failing of this 
model is obvious. The residuals near (—3, +16) 
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Fig. 5. — Residual velocities after removal of 
a model Hubble flow. Residual velocities which 
are less than 1— a are shown as gray; greater than 
l—o" are shown as black (approaching) or white 
(receding). For clarity only the group residual is 
shown for galaxies in groups. 

are the Virgo "s-wave", positive peculiar veloci- 
ties on the near side and negative on the far side 
of Virgo. 

Another way of looking at the quality of the 
fit is to plot observed peculiar velocity as a func- 
tion of model peculiar velocity. In order not to 
be overwhelmed by galaxies which have a large 
virial velocity (which does not appear in the 
mean model velocity) , we plot only those galaxies 
where the model velocity uncertainty is less than 
250 km s -1 . There appears to be a correlation in 
Figure ||] between model and observed, but the 
agreement is not good. 

4.3. Virgo Attractor Infall 

Adding an extended "Virgo Attractor" (VA) 
with 7 = 2 near the location of the Virgo cluster 
at (—3.1, 16.6, —1.6) Mpc, we find an infall ve- 
locity at the Local Group of uy\ = 101 km s -1 , 
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Fig. 6. — Observed peculiar velocities are com- 
pared to model peculiar velocities. Only galax- 
ies with model velocity uncertainty less than 
250 kms -1 (essentially [(H 5r) 2 + 5v 2 1T - m] \ 1 / 2 ) are 
shown, thereby excluding galaxies where a signif- 
icant virial velocity is expected and showing the 
quality of match between the model and observed 
large scale flow field. 

a Virgo thermal velocity dispersion of ova = 
676 km s -1 , and a cutoff radius ryA = 9.3 Mpc. 
The Hubble constant rises slightly to Hq = 72, 
and w = (-320, 195, -82) in the CMB frame. 
This position for the VA is the best fit location of 
a somewhat more elaborate model which also in- 
cludes a Great Attractor. By itself the Virgo At- 
tractor settles to a location of (—3.7, 15.4, —1.7). 
The likelihood is J\f = 345.0, which is a signif- 
icant improvement for just three additional pa- 
rameters. Figure [7| shows the residuals from this 
model. The "s-wave" of Figures || and |H| has been 
removed; the black and white points are evenly 
distributed around the center of the model at- 
tractor at (—3, +16). However, a quadrupole sig- 
nature remains evident, by which residuals in the 
±SGY direction tend to be blueshifted and resid- 
uals in the ±SGX direction tend to be redshifted. 



15 




-40 -20 20 



SGX 

Fig. 7. — Residual velocities after removal of a 
model consisting of Hubble flow, Virgo Attractor, 
and constant velocity vector. 

Figure || shows the galaxy by galaxy residuals 
from this model. While they are somewhat bet- 
ter than Figure |(| the residual quadrupole causes 
substantial scatter. 

The x% = 1-26 value for this fit demonstrates 
that the model is a better match of the veloc- 
ity field than a pure Hubble flow, but the evi- 
dent quadrupole and the peculiar velocity in the 
CMB frame which misses pointing at the Cen- 
taur us cluster by only 10 degrees, indicates a 
new flow in that direction. Following Lilje et al. 
(1986) we can add a quadrupole component to 
this fit, centered on the origin and with no cutoff. 
We find that J\f improves substantially to 320.4, 
and the rest of the parameters change slightly to 
H = 75, n V A = 132, a Y A = 649, r VA = 14.4, 
and w = (-265,200,-97). The quadrupole 
has a stretch of 6.5 km s" 1 Mpc -1 in roughly the 
+X/— Z direction, a stretch of 4.1 kms -1 Mpc -1 
in the +X/+Z direction, and a compression of 
— 10.6 km s -1 Mpc -1 in the ±SGY direction. 
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Fig. 8. — Comparison between model and ob- 
served residual velocities consisting of Hubble 
flow, Virgo Attractor, and constant velocity vec- 
tor. 

4.4. Virgo and Great Attractor Infalls 

The above suggests that it might be profitable 
to try a fit with a second attractor. This is 
our "Great Attractor" (GA) component, and we 
will allow it to vary in both amplitude and po- 
sition. We fit tU+VA+GA in the CMB frame 
and find Hq = 73.5, VA parameters of uva = 
127, <tva = 667, rvA = 14.2, GA parame- 
ters of uqa = 199 and tqa = 28.3, centered 
on (-37.7, 13.3, -18.2) ± (1.7, 2.6, 1.6) Mpc, and 
w = (-146,143,-15) ± (79,34,51) kms" 1 . The 
likelihood has improved to N = 294.9 for 273 
DOF. Figure |9| shows the residuals when this 
model is removed from the observed velocities, 
and figure |l(] shows the residual comparison. 
Like the previous one, this model uses power law 
exponents of 7 = 2 and core radii r c = 2. 

This two attractor plus dipole model appears 
to be a good description of the observed velocity 
field. It constitutes a significant improvement 
over a simple Virgo infall model. If we force the 
overall dipole to zero in the CMB frame (w = 0), 
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Fig. 9. — Residual velocities after removal of a 
model consisting of Hubble flow, Virgo and Great 
Attractors, and constant velocity vector. 

we find that the cutoff radii want to go to infinity, 
and that M rises to 305.3 for 278 DOF. The other 
parameters change quite a bit: Hq = 79.9, uy A = 
207, ova = 675, tva = °o, uq A = 397 and r G A = 
oo, centered on (—37.0,16.8,-16.7) Mpc. This 
illustrates how sensitive some of the parameters 
(such as GA infall velocity) are to an assumed 
velocity reference frame. 

4.5. VA, GA, and a Quadrupole 

Finally, we return to the notion that influences 
from outside the volume where we have data will 
manifest themselves to lowest order dipole 
velocity and a quadrupole, and that our spher- 
ical flow fields are overly simplistic. We see a 
persistent w y velocity and Figure || still seems to 
have a dipole pattern. We now fit a model which 
includes a velocity dipole, the Virgo and Great 
Attractors with free amplitudes and r cut , and a 5 
component quadrupole which is centered on the 
origin and has a cutoff radius of Tq Ua( ^ = 50 Mpc. 
We use a power law exponent of —1.5 for the 
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Fig. 10. — Comparison of velocity residuals af- 
ter removal of a model consisting of Hubble flow, 
Virgo Attractor, GA, and constant velocity vec- 
tor. 

VA and —2.0 for the GA, as suggested by the 
confidence contours illustrated below, and Qm = 
0.2. The positions of the VA and GA are fixed 
at (-3.1,16.6,-1.6) and (-36.7,14.6,-17.1) re- 
spectively (their best-fit locations). We find 
that H = 78 ± 3, u yA = 139 ± 48, r yA = 
12 ± 6, u GA = 289 ± 137, r GA = 50 ± 44, 
and w = (-55,143,-8) ± (102, 41, 62) kms" 1 . 
The likelihood has improved to N = 269.2 for 
272 DOF. Figure |ll] shows the residuals when 
this model is removed from the observed veloci- 
ties, and Figure |l2] shows the residual compari- 
son. The quadrupole found by this model 
has an expansion of 14.6 kms -1 Mpc -1 in the 
(—0.56, —0.16, +0.82) direction, a compression of 
-11.8 km s^ 1 Mpc^ 1 in the (-0.06, +0.99, +0.15) 
direction, and a small compression of -2.8 km s- 1 Mpc 
in the (+0.83, -0.04, +0.56) direction. It and the 
peculiar dipole w y = +143kms _1 are significant 
at more than 3-sigma. 

The reason that introducing a traceless quadru- 
pole changes the model Hq is that the sam- 
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Fig. 11. — Residual velocities after removal of 
a model consisting of Hubble flow, Virgo and 
Great Attractors, constant velocity vector, and 
quadrupole. 
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Fig. 12. — Comparison of velocity residuals af- 
ter removal of a model consisting of Hubble flow, 
Virgo and Great Attractors, constant velocity 
vector, and quadrupole. 



pie galaxies are not isotropically distributed. In 
particular, the major compression axis of the 
quadrupole is in the direction of Virgo, which 
is obviously rich in galaxies. Thus, without the 
quadrupole component, the model is forced to 
compensate for this extra compression by low- 
ering its value for Hq. It is possible that the 
quadrupole arises locally because discrete, spher- 
ical attractors are being used to model flows in 
what is actually an anisotropic potential. We 
discuss this idea further in §[|^. Alternatively 
the dipole and quadrupole could come from un- 
modelled attractors outside our volume. We in- 
vestigate below whether some of the well-known 
galaxy superclusters such as Perseus, Coma, and 
Shapley could be the origin, but we will defer 
a more complete investigation for a subsequent 
paper. 



4.6. Virial Motions 

If we allow the cosmic thermal velocity disper- 
sion to vary from the fixed value of 187 kms -1 
that we have been assuming, then it now set- 
tles in at dcosmic = 187 ± 14 kms -1 for our 
to+VA+GA+Q model. The model likelihood 
does not change in any significant way, of course. 
This result is completely insensitive to the mo- 
tions of galaxies within the Virgo cluster, as these 
are handled by the fitted cluster virial disper- 
sion. Actually, this best-fit "thermal" disper- 
sion also includes a component from the veloc- 
ity measurement errors, but these are relatively 
small and correcting for them puts <r coam i c around 
180 km s" 1 . 

Our fitted value for the Virgo velocity disper- 
sion of <7va ~ 650 km s -1 is identical to those 
found in kinematical studies. The agreement is 
expected, since the velocity data are largely the 
same. For instance, Girardi et al. (1996) con- 
cluded a = 640^g5kms~ 1 for galaxies within 
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2 Mpc of the Virgo core. If we delete the 12 
galaxies within 2 Mpc of the center of the Virgo 
cluster, the fit parameters for the (7va = 1-5) 
model become Hq = 78, uva = 139, -uga = 290 
and w = (—53,142,-8). This is virtually iden- 
tical to the model that includes the Virgo core, 
and demonstrates that our models are able to ac- 
cept large virial velocities and still reliably follow 
large-scale flows. 

4.7. Uncertainties, Covariances, and Con- 
straints 

We have deliberately tried to provide a great 
deal of flexibility in the modelling (for example 
by always fitting for a velocity reference frame). 
However, there are so many parameters making 
up this model that it can be deceptive to quote 
a formal error on a parameter without disclos- 
ing the covariances it has with other parameters. 
The greatest covariances come from the r cut pa- 
rameters, interacting with 7 and infall amplitude 
especially. The next most troublesome covari- 
ances come from the SGX component of the ref- 
erence frame velocity vector w mixing with the 
GA infall amplitude and Hq, and there is also 
a little bit of direct covariance between Hq and 
uqa- Apart from these, the formal covariances 
between parameters are less than 0.6 in magni- 
tude. 

The contour diagrams of this section illustrate 
these covariances as well as graphically show- 
ing the uncertainties in other interesting param- 
eters. The shaded regions of these plots enclose 
the 68, 90, and 95 percent confidence regions for 
two variables, and the bars on the axes show 
the ±1(7 limits for each variable considered sepa- 



rately. Figure |13| shows the joint confidence con- 
tours on the distance, infall velocities, and lo- 
cation on the sky of the Virgo Attractor, and 
Figure |lj| shows the same thing for the Great 
Attractor. The locations of these attrac- 

tors are evidently quite well constrained purely 
by the velocity data. The best-fit distance for the 
Virgo Attractor of 17.1 Mpc is coincident with 
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Fig. 13. — Joint confidence contours of the fitted 
infall velocity, distance, and location of the sky 
of the Virgo Attractor. The points show the lo- 
cation of the galaxies in the RC3 with Vh < 2800. 
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Fig. 14. — Joint confidence contours of the fit- 
ted infall velocity, distance, and location of the 
sky of the GA attractor. The points are from 
the SPS survey; Centaurus lies at (156,-12), 
Abell 1060 is the cluster at (139,-37), and the 
Galactic plane is evident on the left side. 

both the location of NGC 4486 and the median 
of the galaxies in the core of the Virgo cluster 
to within the errors. The location of the GA is 
similar in distance to the most distant galaxies 
attributed to the Cen-30 cluster but lies about 



15 degrees from Cen-30. As figure 14 illustrates, 
this is a complicated region. 

Figure ^ shows how the SGX component w x 
of the dipole velocity affects both the GA in- 
fall amplitude and the Hubble constant. The 
cosmic thermal velocity is little affected by any 
other parameter, and is quite well constrained at 
187 km s . The Hubble constant Hq has a for- 
mal covariance of 0.6 with the GA infall ampli- 
tude, but such a small covariance scarcely affects 
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Fig. 15. — Joint confidence contours of the SGX 
component w x of the dipole velocity with GA 
infall amplitude (left) and Hq (right). 
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Fig. 16. — Joint confidence contours of the cutoff 
radius r cut and power law slope 7 (left) and infall 
amplitude (right) for the Virgo Attractor. 

the determination of each independently. 



Figures 16 and 17 show the confidence con- 
tours for the r cut parameter on the power law ex- 
ponents and the infall velocities of the Virgo and 
Great Attractors. Evidently these parame- 
ters are very poorly constrained individually and 
values of < 7 < 2 for VA and 1 < 7 < 3 for 
GA can be offset by appropriate changes in r cut 
to yield acceptable models. Nevertheless, overall 
run of peculiar velocity with position varies little 
among these models. 



Figure 18 shows the run of net model velocity 
as a function of position along vectors towards 
the Virgo and Great Attractors. The four curves 
are 7 = 0.1,1,2,3, and r cut = 5,9, 19, 00 (VA) 
and 7 = 1,2, 2.5, 3, and r cut = 15, 50, 00, 00 (GA) 
for which J\f is virtually constant around 269. 
(The enclosed density profile of the Virgo At- 
tractor can be described quite accurately as an 
exponential with scale length 5 Mpc.) Except in 
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Fig. 17. — Joint confidence contours of the cutoff 
radius r cut and power law slope 7 (left) and infall 
amplitude (right) for the GA attractor. 

the cores of the attractors, the curves differ very 
little. We also plot the galaxies which lie within 
25° of the vectors (| cos 9\ > 0.9) to give a sense of 
the number of constraining points we have near 
the attractors (although the attractor models are 
affected by galaxies all over the sky) . The agree- 
ment between the galaxy and model velocities in 
this plot should only be very approximate close 
to the attractors because the angle between the 
line of sight and the vector from the galaxy to 
the attractor becomes much greater than 25°. 



Figure 19 illustrates this slightly differently. 
These panels show how "Ui n f a ii for each of the two 
attractors individually varies as a function of dis- 
tance from the attractor, given these four mod- 
els with very different 7 's. The models diverge 
close to and far from the attractors but there is 
a broad range where the models give substan- 
tially the same values for infall velocity (6 < r < 
25 Mpc for the VA and 10 < r < 50 Mpc for 
the GA). We draw power-law approximations for 
these two attractors on the plots whose slopes 
are —1.1, although it is apparent that all mod- 
els for both attractors curve away from a pure 
power law over these radial ranges. The ratio of 
attractor infall velocity to Hubble flow is for the 
Virgo Attractor 
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Fig. 18. — The net model velocity is shown 
as a function of position along vectors running 
through the Virgo and Great Attractors using 
various values of 7 and r cut which leave J\f 
constant. The dashed line is just H^r. The 
points are those galaxies within 25° of the vectors 
(|cos0| > 0.9). 



and for the Great Attractor 

( r 



0.09 



V^GA 



-2.1 



(14) 



where ryA and tqa are the distances of the Local 
Group from the attractors. 

Not surprisingly, for a given value of Qm we 
find quite consistent values for 5, the overden- 
sity within our radius, for these models. For the 
Virgo Attractor, values of 0.1 < 7 < 3 lead to 
values of 6 = 1.0 ± 0.07, or M VA = 7 x 10 14 M . 
However, if the quadrupole is caused by a non- 
spherical Virgo (as we suggest below), these es- 
timates for 5 and Mva are not meaningful, al- 
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Fig. 19. — The model infall velocity is shown as 
a function of distance from the Virgo and Great 
Attractors for the four sets of 7 values. A solid 
line delineates the range where there is substan- 
tial agreement and where we believe the models 
are well constrained by the data. 



though Equation (13) is. For the Great At- 
tractor, models with 1 < 7 < 3 lead to values 
of 5 = 0.73 ± 0.14, or M GA = 8 x 10 15 M . 
(These masses, of course, are the masses in ex- 
cess of background density within spheres cen- 
tered on the attractors with radius reaching the 
Local Group.) From this it is easy to solve for 
the radius and angle subtended by the zero veloc- 
ity surface (sphere where infall cancels outward 
Hubble flow) around these attractors. For the 
Virgo Attractor this is 6 Mpc, or 20°; for the GA 
it is 14 Mpc, or again about 18°. We see from 
Figure that by 5 « 1 (u/H r xi 0.1) the mod- 
els do not fit well unless they are falling at least 
as steeply as r -3 . In particular, pure isother- 
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Fig. 20. — Contours are shown of our model of 
the radial component of the flow field in a plane 
which cuts through the Local Group, Virgo, and 
Great Attractors. The model velocity of the 
Local Group with respect to the CMB frame 
is apparent in the discontinuity at the origin. 
The Virgo Attractor is found at "SGX,SGY" = 
(—3, +17) (quotes because the plane doesn't cor- 
respond perfectly to supergalactic coordinates); 
the Great Attractor is at (—42, +16), and radial 
lines are drawn at ±10° and ±20° from it. 

mal models (7 = 2, unmitigated by any cutoff) 
would provide poor representations of the mass 
distributions. 

The GA distance is well constrained in these 
models at cIqa = 43+3 Mpc. This is a bit behind 
most of the galaxies comprising the Cen-30 and 
Cen-45 clusters. In fact, the Cen-30 "cluster" 
appears to be galaxies distributed over 5-10 Mpc 
on the near side of the GA and the "s-wave" dis- 
tortion of the Hubble flow causes them to share 
similar redshifts. This is illustrated in Figures 20 
and ^l], which show contours of our flow model 
with and without survey galaxies overlaid. We 
do not have a good constraint on the GA ampli- 
tude closer than 10 Mpc, so the very dramatic 



Fig. 21. — This is the same figure as the previ- 
ous one except that our survey galaxies are over- 
plotted. Note that the Centaurus galaxies which 
appear to be going through the Great Attractor 
actually pass above it by about 15°, hence lie in 
the stall zone. 

peculiar velocities inside that radius may well be 
overestimates. In this projection, the Centaurus 
galaxies fall above the GA, but they lie in the 
stall region near ±15° angular separation, and 
the model and data agree on a typical CMB ve- 
locity of 3150 km s -1 . Similarly, the Ursa Major 
group lies in the stall region near the VA with 
a CMB velocity of 1150kms _1 , and the Virgo 
Southern extension lies in a stall zone with CMB 
velocity of 1500 km s -1 . The general trend of 
the model for i>cmb to range from 1100 kms" 1 
to 1500 kms" 1 as SGL swings from 85° to 125° 
is apparent in redshifts surveys of the region, as 



is illustrated in Figure |22, 

We note that the GA distance increases by 
about 10% when the 3 high velocity Cen-45 
galaxies (with good quality data) are included 
in the modeling. However, the likelihood of the 
model suffers badly in trying to accomodate just 
these 3 additional DOF, rising from N = 269 
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Fig. 22. — The CMB velocities of galaxies from 
the RC3 are plotted function of angular sepa- 
ration from the GA (left), and function of su- 
pergalactic longitude SGL (right). The concen- 
tration at f cmb = 3150 (left) and the trend for 
vcmb to change from 1100 to 1500 km s -1 (right) 
inherent in our model is apparent in these red- 
shift data. 

to N = 295. Thus, for now we have cho- 
sen to exclude them. A detailed discussion of 
the distribution and motions of galaxies in the 
GA/Centaurus region, and a comparison of the 
SBF flow model results with those of previous 
GA models will be given by Dressier et al. (1999). 

4.8. Influence of Perseus-Pisces and Shap- 
ley 

We now ask if the Perseus-Pisces superclus- 
ter and/or the Shapley Concentration could be 
sources of the peculiar dipole w and quadrupole 
that we find. To test this idea, we add to our 
standard tu+VA+GA+Q model compact attrac- 
tors representing Perseus-Pisces at (65, —15, —17) 
and Shapley at (-170,120,-30), allowing the 
masses to be free parameters. The likelihood M 
improves to 267.2, which is not a significant im- 
provement for two new parameters over the pre- 
vious value of 269.2. For Perseus-Pisces, we find 
an insignificant velocity amplitude of v amp — 

59± 

314 km s ; recall that this is normalized at the 
fiducial radius of 50 Mpc. For Shapley, we find 
an extraordinary u amp = 17926 km s , which for 
the assigned distance of 210 Mpc corresponds to 
a velocity of 1013kms _1 at the Local Group. 
The other parameters change a great deal to 
offset this, especially w which adjusts itself to 



(920, -530,210). Equally telling however, is the 
fact that the quadrupole grows to have eigenvec- 
tors of 22, —20, and —2 kms _1 Mpc _1 , indicat- 
ing that the inclusion of Perseus and Shapley has 
made this model less plausible. 

There is a great deal of covariance between the 
distant Shapley attractor and the w term, so we 
ran another model including the Shapley attrac- 
tor but with w set to zero. This model gave a ve- 
locity amplitude w amp = 2533±717kms~ 1 , corre- 
sponding at the Local Group to 144 ± 41 km s -1 . 
The likelihood of this model is N = 274.0, which 
for only one additional parameter is a substan- 
tial improvement over the N of 278.5 found for 
a model with VA+GA+Q but no w or Shapley 
component. However, the fit is worse than our 
standard model with w and no distant attractors, 
although the difference in the fitted quadrupole 
components for these two models is now insignif- 
icant. (Adding a Perseus attractor gives no fur- 
ther improvement; it is in the wrong direction to 
compensate for the lost flexibility of w.) 

We conclude that Perseus-Pisces exerts no sig- 
nificant influence on our local volume, consis- 
tent with previous studies in this direction (e.g., 
Willick 1990; Courteau et al. 1993; Hudson et al. 
1997), and that the Shapley Concentration is not 
the source of the peculiar quadrupole in our mod- 
eling, although it might contribute some part of 
the observed dipole. Bearing this last point in 
mind, we continue with our standard model be- 
cause it provides a better fit to the data than the 
model which replaces w with a Shapley attrac- 



tor. In addition, we show below in §5.2 that the 



dipole can be removed by simply translating the 
quadrupole origin to the center of Virgo. 

Similarly we experimented with including nearby 
Abell clusters, selected under the assumption 
that their galaxy counts are proportional to their 
masses and ranking them by M/r 2 . Using a suite 
of the 22 most significant ones (including Coma, 
Perseus, A2199, etc.) and allowing their mass- 
to-luminosity ratio to be a free parameter gives a 
likelihood N = 268.8, not a significant improve- 
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ment. The model parameters are virtually un- 
changed, including the quadrupole and the dipole 
velocity. 

4.9. n M 

Lacking any a priori information about the 
values of the 6's, we have very little constraint on 
Only the non-linear component to Yahil's 
p 1 / 4 law (Eq. provides any leverage at all. Nev- 
ertheless, if we permit VLm to be a free parameter 
the best fitting model comes in with N = 269.2, 
and VLm = 0.08 ± 0.33. If we refit our standard 
model with Qm = 1, we find N = 269.5, again in- 
dicating how little direct handle we have on £Im- 
The other parameters change very little, but as 
expected, the combination of S Oli^ 3 for the two 
attractors is relatively constant for the various 
values of VLm at 0.33 for the Virgo Attractor and 
0.27 for the Great Attractor. We plan in a fu- 
ture work to compare our measured peculiar ve- 
locities to the galaxy distribution and so obtain 
constraints on £Im- 



5. Discussion 

5.1. Mass Distributions of the Virgo and 
Great Attractors 

We find good fits to the velocity field us- 
ing attractors whose density distributions go as 
p ~ r~ 2 when r is small enough that 5 > 1, 
but steepening to p ~ r -3 by S « 1. Most previ- 
ous parametric flow models have adopted isother- 
mal 7 = 2 attractors (e.g., Lynden-Bell et al. 
1988; Han & Mould 1990), as 7 was not well- 
constrained by their data sets. Aaronson et al. 
(1989) did allow 7 to vary in their model, al- 
though they constrained it to be the same for 
both the Virgo and Great Attractor mass con- 
centrations, and found 7 ~ 2. Faber & Burstein 
(1988) also experimented with different 7's and 
found that Virgo was roughly isothermal, and the 
more distant GA preferred 7 ~ 2.7. This is sim- 
ilar to what we find with r cut set to infinity and 
7 allowed to vary, but we prefer the approach 
adopted in the previous section because of the 
substantial improvement of the fit. 

The isothermal cluster assumption was moti- 
vated by dynamical evidence that the enclosed 
mass in clusters increased approximately linearly 
(e.g., Press h Davis 1982). Observations con- 
tinue to show flat velocity dispersion profiles for 
many clusters, in support of this assumption 
(e.g., Fadda et al. 1996). Weak lensing obser- 
vations are also consistent with isothermal, or 
somewhat steeper, mass distributions (Squires 
et al. 1996a, 1996b). However, none of these 
studies probe beyond a few Mpc from the clus- 
ter centers, whereas we are dealing here with the 
distributions over tens of Mpc. 

Numerical studies by Navarro, Frenk, &; White 
(1996, 1997) suggest that cluster dark matter 
halos in hierarchical collapse models should go 
as r -3 at large radii. Earlier numerical work 
on the collapse of isolated halos produced even 
steeper r -4 profiles at large radii (Dubinski & 
Carlberg 1991). The total mass of course de- 
pends strongly on 7, with the 7 < 3 yielding an 
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infinite mass; thus at some point the distribution 
must steepen. Our observed mass profiles there- 
fore are quite reasonable in terms of flat velocity 
dispersion profiles and the numerical simulations. 

5.2. Dipoles and Quadrupoles 

Although we find a 4-er, 150 kms -1 dipole 
residual pointing in the SGY direction, we re- 
gard this result as an indication that our sur- 
vey volume, which extends to about 3000 km s -1 
or 40Mpc, is essentially at rest with respect to 
the CMB. Indeed, we can clearly create or re- 
move any dipole velocity we like by shifting the 
origin of the quadrupole contribution, as is evi- 
dent from Equation (||), provided the quadrupole 
matrix is invertible. Applying the inverse of 
the quadrupole to the dipole velocity w of our 
best fit model implies a quadrupole origin of 
(—17,13,-9), rather close to the VA, although 
moved somewhat towards the GA. If we shift the 
quadrupole origin from the Local Group to the 
center of the Virgo Attractor (which can then 
partially mimic the effects of a flattened attrac- 
tor), choose a cutoff r qua d = 35 Mpc, and fit with 
w = 0, we find N = 272.0, which is only 3 greater 
with effectively three fewer degrees of freedom; it 
is 6.5 less (with the same DOF) than the model 
with w = and the quadrupole centered on 
the LG. Alternatively we can fit for w with the 
quadrupole centered on the Virgo Attractor, and 
we find w = (—91, —7, —63) or within 1-cr of zero 
in each coordinate. We therefore think it likely 
that the galaxies in our survey volume are on av- 
erage at rest with respect to the CMB, and that 
the dipole and quadrupole terms of our standard 
model are acting as lowest order correction terms 
to our spherical Virgo Attractor. 

We found previously that addition of a quadrupole 
component to the two attractor model caused 
the Hubble constant to increase from 73 to 78 
as the likelihood improved. This is illustrated 
in the left panel of Figure ^ which shows confi- 
dence contours of our standard model as a func- 
tion of Hq and quadrupole amplitude relative to 
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Fig. 23. — Joint confidence contours of the Hub- 
ble constant and relative quadrupole amplitude 
when the quadrupole is centered on the Local 
Group are shown on the left. Confidence con- 
tours of Hubble constant and quadrupole cutoff 
radius when the quadrupole is centered on Virgo 
instead are displayed on the right. 

the best fit value. We find that moving the ori- 
gin of the quadrupole to Virgo also has a slight 
side effect on the value of the Hubble constant. 
If there were no cutoff in the quadrupole its ori- 
gin would be completely degenerate with w, but 
when when the quadrupole is centered on Virgo 
we find that Hq drops slightly to 76, and there 
is some covariance between Hq and 7" qua <} in the 
sense that larger r qua d leads to larger Hq, as seen 
in the right panel of Figure 23. Because this 
use of a quadrupole correction to a spherical in- 
fall model seems like a poor approximation to 
a flattened potential, and because the eigenvec- 
tors of the quadrupole are tipped by 45° from 
the SGZ symmetry axis of the galaxy distribu- 
tion, we do not prefer this model over the one 
which has the quadrupole centered on the Local 
Group. Nevertheless, we are mindful of the sen- 
sitivity of Hq to this possibility and we therefore 
choose Hq = 77 ± 4 as the most likely range we 
can derive from these data and models. 

Unlike Lynden-Bell et al. (1988), we do not 
find that Centaurus has a much larger radial pe- 
culiar velocity in the CMB frame than the Local 
Group. Instead it appears to have a significant 
extent falling into the attractor at 43 Mpc and at 
least some galaxies (e.g., NGC 4767) have no ra- 
dial velocity with respect to the CMB. Dressier 
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et al. will discuss this issue in greater depth, 
but it appears that the flow in the Great Attrac- 
ted region is fairly complex and may need many 
more high accuracy distances before it can be 
completely untangled. 

5.3. Motion of the Local Group 

The Local Group is moving at 627 kms" 1 , 
with components (—406, 352, —324) kms -1 , with 
respect to the CMB frame, but we have not used 
this fact at all in our modelling. To what ex- 
tent can we account for this motion of the Lo- 
cal Group reference frame in the context of our 
model? The three contributors to the model ve- 
locity at the Local Group are as follows: 
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We can account for w x and w y quite well; the 
residual of —79 kms -1 is not statistically signif- 
icant and is in any case a reasonable peculiar 
velocity for the Local Group with respect to the 
nearby Hubble flow. 

The w z residual of — 188 kms -1 is much too 
large to be an unsuspected motion of the Milky 
Way with respect to the LG frame. The models 
by Yahil et al. (1977) and Courteau & van den 
Bergh (1999), among others, require that this be 
a motion of the entire Local Group. Our sample 
does not have any nearby galaxies which are far 
enough out of the supergalactic plane to ascer- 
tain whether this motion is restricted to the Lo- 
cal Group or extends farther. The nearest galaxy 
in our sample with |SGB| > 20° is 8 Mpc distant, 
and at SGB= —33° we could not distinguish be- 
tween v z = and v z = —188. We note that this 
velocity is similar in magnitude to our best-fit 
1-d thermal velocity. 

This peculiar motion of the Local Group in 
the general direction of — SGZ has been noted be- 
fore and was named the "Local Anomaly" (LA) 



by Tully (1988) and Faber & Burstein (1988), 
who both attributed it to inhomogeneities in the 
galaxy distribution within cz ~ 1000 kms" 1 , in- 
cluding the "Local Void" (Tully k Fisher 1987). 
These authors found fairly large LA velocities: 
the Faber Sz Burstein value was 360 kms" 1 to- 
wards (/, b) = (199°, 0°), or (169, 30, -317) kms" 1 
in supergalactic components. Han & Mould 
(1990) modeled the LA in more detail and found 
a smaller net motion of 236 kms" 1 towards (/, b) = 
(205°, 11°), or (88, 68, -208) kms" 1 . All three of 
these sets of authors used the Aaronson et al. 
(1982b) infrared Tully-Fisher data set to measure 
the LA. Aaronson et al. (1982a) themselves found 
a LG peculiar motion of (38,86,-150) kms" 1 
and attributed the large —SGZ velocity to a dis- 
placement of the Local Group above the super- 
galactic plane and a resultant downward acceler- 
ation. 

Lahav et al. (1993) suggested that the LA 
could be explained by the combined effects of 
the Local Void in the +SGZ direction and the 
Puppis concentration of galaxies towards —SGZ, 
with possible smaller contributions from Fornax 
and Eridanus. The Puppis concentration was re- 
vealed by Hi surveys of IRA S-selected galaxies; 
it is a rich, but loose, association of galaxies lying 
in the Galactic plane at (l,b) = (245°, 0°) with 
a mean velocity of about 2000 kms" 1 (Yamada 
et al. 1994). Lu & Freudling (1995) surveyed the 
region behind the Galactic plane towards the di- 
rection of the Faber & Burstein (1988) and Han 
& Mould (1990) LA vectors and disputed the 
conclusion of Lahav et al. Since they found no 
other significant excesses of nearby galaxies, Lu 
& Freudling concluded that any motion arising 
from a gravitional pull of nearby galaxies and a 
push from the Local Void, whose center is nearly 
diametrically opposite Puppis, would have to be 
directed almost straight towards Puppis, nearly 
45° from the supposed direction of the LA vector. 

In fact, the Local Anomaly we find here has 
an apex at (/, b) = (250°, —13°), remarkably close 
to the direction of the Puppis concentration and 
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to the anticenter of the Local Void as sketched 
out by Lu & Freudling. However, the magni- 
tude we find is significantly smaller than the 
Faber & Burstein value adopted by Lahav et al. 
(1993); thus, the ~ 50 kms -1 contribution they 
estimated for the Puppis galaxies becomes less 
necessary. Likewise, Han & Mould (1990) con- 
cluded that the Local Void might account for 
their fairly modest LA motion, if the mean posi- 
tion of the Void was a free parameter. 

We experimented with adding a Gaussian void, 
as defined by Eq. (||), to our standard model 
along the positive SGZ axis. We first looked at a 
very small Local Void with both the distance and 
Gaussian size set to 10 Mpc and fitted for the 
velocity amplitude. Since there is some covari- 
ance between f am p and w z , we set w z = 0. We 
find a best-fit amplitude of u amp = —8.5 ± 5.6 
kms -1 , and an improvement in J\f from 269.2 
to 267.7, marginally significant. This void pro- 
vides a push of —65 kms -1 at the location of the 
LG. The model parameters show little change, 
although Hq drops from 78 to 77 kms -1 Mpc -1 . 
Alternatively, we tried adding a Local Void us- 
ing the distance and size of 20 Mpc suggested by 
redshift surveys and fixing the amplitude at the 
value of Vamp = —88 kms -1 which provides the 
push of —188 kms -1 at the Local Group that we 
are trying to understand. The central 5 of this 
void is close to —1. This time we do allow w z 
vary. The fit improves still more to N = 265.9 
and w z is still negligible at 20 kms -1 , but the 
amplitudes of the VA and GA shift quite a bit 
to 172 and 159 kms -1 respectively, w x and w y 
also move to —136 and +151 kms" 1 , and Hq 
decreases to 73 kms -1 Mpc -1 . Obviously fine- 
tuning the free parameters could produce an even 
better match to the excess LG motion, but this 
is precisely the sort of ad hoc modelling which 
we think is not appropriate. This void's proper- 
ties are not well-constrained in our model nor are 
there well- motivated external constraints, so we 
choose not to make this adjustment to our stan- 
dard model, but simply note that our peculiar 



velocity data do seem to support the idea that 
there is a 20 Mpc void above the Local Group 
which is nearly empty at its center. 

Another interesting issue is the peculiar ve- 
locity of the Local Group with respect to the 
Virgo Attractor, since this was the historical 
means of determining the Virgo infall amplitude. 
The Local Group is approaching Virgo at about 
139 kms -1 due to the Virgo Attractor itself; the 
Great Attractor causes the LG and Virgo to con- 
verge at 97 kms -1 ; the extra quadrupole com- 
ponent results in the LG and Virgo approaching 
each other at an additional 173 kms -1 . The total 
model peculiar velocity between the Local Group 
and Virgo is therefore 409 kms -1 . 

It is unclear how much "infall velocity" to as- 
cribe to the Virgo Attractor. Using our model of 
the Virgo Attractor spherical distribution, 
it only tugs on us at 139 kms -1 . However, if we 
regard the entire quadrupole clS cl component of 
an anisotropic Virgo Attractor, the net tug on 
the LG is 303 kms -1 . These considerations may 
partially account for the widely disparate values 
which have been reported for the Virgo infall am- 
plitude. Depending on the sample one is using 
and whether one's fit is dominated by model am- 
plitudes near our location or averaged all around 
the Virgo supercluster, it is possible to get naive 
"infall velocities" (i.e., peculiar velocity between 
Virgo and the LG) anywhere between 140 kms -1 
and 400 kms -1 . 

5.4. Comparison with Other Flow Stud- 
ies 

Our analysis finds more modest infall veloci- 
ties at the Local Group and a smaller bulk flow 
than most previous studies. Part of the reason 
for our smaller infall velocities is that our mod- 
eling has been more flexible in not requiring a 
single power-law profile for the massive attrac- 
tors. Thus, our model can accommodate fairly 
large infalls close to the attractors without re- 
quiring large motions at the Local Group. In 
addition, we have not required the spherical in- 



27 



falls to account for all of the Local Group mo- 
tion in the CMB frame, but rather have allowed 
the reference frame to be a free parameter. In 
fact, setting w = 0, r cut = oo, and omitting the 
quadrupole term yields a model similar to that 
of Lynden-Bell et al. (1988), but at a significant 
cost in the likelihood (see § [4.4[) . 

Our best-fit bulk flow w of ~ 150 kms -1 to- 
wards (I, b) = (294°, +67°) ± 43°, and the strong 
indications that even this modest value may be 
an artifact of the model, appears inconsistent 
with several recent studies which find large-scale 
streaming motions of amplitude ~ 700 kms -1 for 
volumes encompassing our survey volume (e.g., 
Lauer & Postman 1994; Hudson et al. 1999; 
Willick 1999). However, it is quite consistent 
with the bulk flows found by the /-band Tully- 
Fisher study of Giovanelli et al. (1998) (using 
their "Case b" weighting). These authors report 
bulk flows of 151 ± 120 kms" 1 towards (l,b) = 
(295°, +28°) ±45° for their full sample of 24 clus- 
ters extending out to 9000 kms -1 , and 131 ± 
gOkms" 1 towards (l,b) = (325°, +62°) ± 60° for 
their 17 clusters within 6000 kms -1 . 

Our best-fit 1-d thermal velocity is cr C osmic = 
180 ± 14 km/s, after making allowance for ve- 
locity measurement uncertainty. Lynden-Bell 
et al. (1988) found a best-fit thermal dispersion 
of 250 ± 40 in their favored model and concluded 
that this result was mainly due to small-scale 
flows, not virial dispersions in groups. 

Strauss, Cen, & Ostriker (1993) calculated the 
thermal velocity dispersions from several peculiar 
motions surveys in order to compare the result- 
ing Mach numbers (the ratio of the mean bulk 
flow to the 3-d thermal dispersion) to expecta- 
tions from various cosmologies. Their estimates 
of the 1-d dispersions from the observational sur- 
veys ranged from 250 to 420 kins" 1 . More re- 
cent large-scale flow studies often simply assume 
a value cr COS mic = 250 kms" 1 for clusters and find 
that the resulting x% ~ 1 ( e -S-' Willick 1999; 
Hudson et al. 1999). The Giovanelli et al. (1998) 
cluster study allowed this to be a free parame- 



ter and found cr cosmic — 

300 ±80 kms" 1 . How- 
ever, on small enough scales there is evidence 
that the dispersion of galaxies (field spirals in 
particular) about their local flow field is closer 
to ~125 kms" 1 (Davis, Miller, & White 1997; 
Willick et al. 1997). This value may better re- 
flect the virial velocities of galaxies within small 
groups. 

5.5. The Revised Calibration and Hq 

The revised Mj calibration of Eq. (||) carries 
a formal zero-point uncertainty of 0.10 mag from 
the SBF/Cepheid galaxy comparison. This un- 
certainty is larger than the 0.07 mag quoted in 
SBF-I because that paper used the group com- 
parison, which averaged many SBF distances for 
each Cepheid distance. As stated above, the 
group comparison is 0.13-0.19 mag fainter than 
the galaxy comparison. By choosing the galaxy 
calibration over the group one, we have implicitly 
assumed that the difference is due to a system- 
atic offset in the mean distances of the HST Key 
Project Cepheid galaxies and the groups in which 
they are thought to reside. The alternative ex- 
planation is that the measured SBF magnitudes 
are systematically too bright for the spiral bulges, 
either because of population effects or because of 
extra variance from dust or other features. Fer- 
rarese et al. (1999) have likewise chosen the for- 
mer explanation and have derived a calibration 
similar to that given here. 

With this new calibration and depending on 
the flow model, we have derived values of Hq 
in the range of 70-80 km s -1 Mpc -1 . This range 
correlates with the amount of mass in the model 
attractors. The value of Hq = 70 corresponds 
to a pure Hubble flow plus overall dipole model; 
at the other extreme with Hq = 80 is a dual 
isothermal (7 = 2) attractor model with r cut set 
to infinity. The former model is ruled out by 
the extremely poor x%i the latter model overes- 
timates the attractor masses, thus overestimat- 
ing Hq. When we allow r cut to vary and pro- 
duce an r -3 (and steeper) large radius falloff 
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to the models, Hq comes out near 74. There 
is, however, a significant degree of anisotropy in 
the data with respect to this model as evidenced 
by the dramatic improvement in likelihood when 
a quadrupole component is introduced into the 
model. As it happens, relieving this anisotropy 
with a quadrupole centered on the Local Group 
causes our value for Hq to rise somewhat to 78, 
whereas a quadrupole centered on Virgo yields 
an Hq of 76. Despite our qualms about the ad 
hoc nature of the latter model, the fact that it 
has a w consistent with zero causes us to prefer 
it slightly for estimating Hq. Thus, for our best 
estimate of a flow-corrected but unbiased value, 
we adopt 

H = 77±4±7kms- 1 Mpc- 1 . (15) 

The first errorbar is the quadrature sum of the 
internal uncertainty for a single flow model and 
an estimate of the likely range from the differ- 
ent models. The second errorbar includes the 
0.10 mag uncertainty in the fitted calibration to 
Cepheids and an estimated 0.17 mag uncertainty 
in the Cepheid zero point itself (see Ferrarese 
et al.); these systematic uncertainties also have 
been combined in quadrature. 

Our value for Hq differs from the value of 
69kms _1 Mpc~ 1 found by the Ferrarese et al. 
(1999) analysis of SBF partially because of our 
choice of zero point and partially because of the 
flow models used. We have done extensive mod- 
eling to estimate and remove the effects of large- 
scale flows from our data. Another difference is 
that we have used the full survey data set in de- 
riving Hq, rather than just four SBF measure- 
ments from HST. Consequently, our best value 
is 11% larger than theirs, although marginally 
consistent within the random uncertainty. 

As with the Key Project results, if the Cepheid 
calibration is changed so that the distance mod- 
ulus of the Large Magellanic Cloud is 0.2 mag 
smaller than the standard 18.5 mag value, then 
Hq would increase by 10%; if the Kennicutt et al. 
(1998) metallicity dependence for the Cepheid 



period-luminosity relation is correct, then Hq 
would decrease by 5%. However, our calibration 
would become quite discordant with the theo- 
retical calibration using Worthey's models if the 
zero point moved by more than about 0.2 mag, 
and those models are largely based on the RR 
Lyrae and parallax distance scales, rather than 
Cepheids. 

6. Summary and Conclusions 

We have shown that the SBF technique is a 
powerful tool in mapping the local and large- 
scale velocity field through accurate distance mea- 
surements to early-type galaxies. Judging by our 
best data, for which the product of seeing and 
distance PD < 1.0, distance measurements with 
an accuracy of ±5% are attainable for any early- 
type galaxies with sufficiently good data. Dis- 
tance measurements of this accuracy are essen- 
tially immune to Malmquist-like biases. 

The parametric model we have developed is 
unique in the way it handles both local pecu- 
liar motions and velocity dispersions, particu- 
larly for galaxy clusters. Although it is im- 
portant to follow this first analysis with non- 
parametric models, we find it remarkable that 
nearly all of the deviation from a smooth Hub- 
ble flow can be accounted for by the addition 
of two spherical, truncated isothermal attractors, 
one centered on the Virgo cluster and the other 
centered just beyond the most distant of the el- 
lipticals associated with the Cen-30 and Cen-45 
clusters. These two attractors appear to center 
on enhancements in galaxy distribution, consis- 
tent with the idea that galaxy light traces mass 
in at least some approximate way. The improve- 
ment in the model with the inclusion of an addi- 
tional dipole and quadrupole is statistically sig- 
nificant, but our analysis indicates that these 
terms are more likely the result of inadequacies in 
the spherical attractor model, rather than mani- 
festations of additional mass fluctuations beyond 
our survey volume. Further measurements, par- 
ticularly beyond the R « 40 Mpc limit of most 
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of our data, will be able to discriminate between 
competing explanations, but it is clear from the 
size and direction of these components that the 
influence of mass overdensities in the Perseus- 
Pisces supercluster or Shapley concentration do 
not play a substantial role in perturbing the ve- 
locity field within 3000 km s -1 . 

Not surprisingly, we find strong covariance be- 
tween pairs of parameters, for example, the am- 
plitude of the GA infall and the strength of the 
dipole component in this direction, or the cutoff 
radius and the power-law exponent of mass dis- 
tribution attributed to the attractors. New data 
are unlikely to reduce some of these covariances, 
but additional measurements on the backside of 
the Great Attractor, for example, will be able 
to constrain much better the contribution from 
more distant concentrations that produce a local 
dipole. 

Our standard flow model resolves most of 
the Local Group's motion with respect to the 
CMB, except for the continued presence of a 
nearly 200 km s" 1 component perpendicular to 
the supergalactic plane. This "Local Anomaly" 
has been seen before in other data sets. The 
amplitude we find is generally smaller than in 
those previous studies, and crude modeling of 
this "anomalous" motion as a "push" from a 
nearby void in the +SGZ direction can account 
for nearly all of it. 

Our data are considerably more accurate than 
those used in the analysis of Lynden-Bell et 
al. (1988). Our best fit model, while qualita- 
tively similar, is significantly different, but this is 
mainly a product of our inclusion of other com- 
ponents, for example, the inclusion of a dipole 
in addition to the two attractors. In particu- 
lar, the thermal velocity of our sample is only 
180 km s^ 1 , significantly lower than the Lynden- 
Bell et al. study. A more detailed comparison 
of our study with previous results will appear in 
Dressier et al. (1999). 

We find zero or at most a very small dipole 
motion in the R < 3000 km s^ 1 volume, which 



is consistent with the study of Giovanelli et al. 
(1998), mildly inconsistent with the larger bulk 
flow reported by Courteau et al. (1994), and 
very inconsistent with the large-amplitude, large- 
scale flows reported by Lauer & Postman (1994), 
Hudson et al. (1999), and Willick (1999) un- 
less the R < 3000 km s -1 is at rest with re- 
spect to the CMB while the thick shell outside 
3000 km s" 1 < R < 10, 000 km s" 1 in rapidly 
moving, a situation we think is very unlikely. 

Finally, we provide a new zero point for the 
SBF relation, based on SBF distances to galax- 
ies that have Cepheid distances determinedby the 
HST Key Project. The Hubble constant that 
results from this calibration, Hq = 77 ± 4 ± 
7kms~ 1 Mpc~ 1 is somewhat higher than results 
of that study, partially because we choose a dif- 
ferent zero point and partially because we have a 
more extensive set of galaxies and a more sophis- 
ticated model. We consider the zero point issue 
to be unresolved and in need of further study. 

Our study shows the power of accurate SBF 
distance measures in advancing our understand- 
ing of the local mass distribution and the dis- 
tortions in the velocity field that result from 
it. Extending the SBF technique, through big- 
ger aperture telescopes with better psf perfor- 
mance, through near-IR SBF measurements, and 
through observations with HST and future space 
telescopes, promises to provide decisive answers 
for the long-asked questions regarding the large- 
scale distribution in the local universe. 

We are grateful to all our friends who have 
helped us collect these data over the years. This 
work was supported primarily by NSF grant 
AST9401519. Additional support for this work 
was provided by NASA through grant number 
GO-06579.02-95A from the Space Telescope Sci- 
ence Institute, which is operated by the Associ- 
ation of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. 
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A. Appendix. The flow model 

Our standard flow model consists of 5 com- 
ponents: Hubble flow, peculiar velocity, Virgo 
Attractor, Great Attractor, and quadrupole. Al- 
ternative models are possible as discussed above 
which center the quadrupole on Virgo or use dif- 
ferent values of 7, but they give very nearly the 
same model velocity as a function of position. 
One caveat is that our measure of the mean 
flow velocities around the VA and GA are ex- 
tremely uncertain for r < 6 Mpc for the VA and 
r < 10 Mpc for the GA. 

We work in supergalactic coordinates, using 
a transformation from galactic to supergalactic 
coordinates by applying the rotation matrix 

(— sin 4>s + cos (f>s \ 

— sin 6s cos 4>s — sin 9s sin 4>s cos 9s 
+ cos 6s cos 4>s + cos 9s sin ips sin 9s J 

(Al) 

to galactic cartesian coordinates, where 4>s = 
47.37° and 9 S = 6.32°. 

Our value for the Sun's motion with respect 
to the CMB (hence transformation from helio- 
centric velocities to CMB reference frame) comes 
from Lineweaver et al. (1996), who find the sun 
moving at 369 km s" 1 towards SGL = 264.31° 
and SGB = 48.05°. 

Given a position (x,y,z) in supergalactic co- 
ordinates with respect to the Local Group, the 
first two terms of our model are 

v = H f+ w (A2) 

where H = 78.4 and w = (-55, 143, -8). 

The two attractors have a radial infall around 
them. Each causes a deviation from the above 
velocity of the form (equation |j) 

Uinfaii = \h t a nlf 5(r A ) {l+SirA))- 1 ^ (A3) 

where r A is the distance from the point of interest 
(x,y,z) to the attractor, Q,m = 0.2, and 8(r A ) is 
the overdensity in units of the background den- 



sity at distance r A from the attractor. The over- 
density 5(r A ) is calculated from equation 0: 

""» = T^r(^) h^ 73 - 1 • 

(A4) 

5o is evaluated by plugging in the distance from 
the LG to the attractor and matching to the 
model overdensities given below. The contribu- 
tion to v is obtained by resolving the radial Ui n f a u 
into components along f—f A . 

The Virgo Attractor is found at (—3.1, 16.6, —1.6) 
Mpc; it has overdensity 5y = 0.974; its infall ex- 
ponent is 7y = 1.5; the core radius is r c = 2; and 
the cutoff radius is r cut = 11.7. 

The Great Attractor is found at (-36.7, 14.6, -17.1) 
Mpc; it has overdensity 5q = 0.781; its infall ex- 
ponent is 7g = 2.0; the core radius is r c = 2; and 
the cutoff radius is r cut = 49.5. 

Finally, the quadrupole contribution is found 
by multiplying r by the matrix 

/ 2.57 2.04 -7.87 \ 
e - r2 /Ku*d 2.04 -11.12 -3.63 , (A5) 
V-7.87 -3.63 8.55/ 

where the quadrupole cutoff radius is r qua d = 
50 Mpc. 

Our model for the (1-d) velocity dispersion 
as a function of position has four components 
added in quadrature. The first is an overall 
thermal velocity of 187 km s -1 . The next two 
are components centered on the VA and GA 
at the positions listed above, with amplitudes 
650 km s -1 and 500 km s -1 respectively at the 
center, but falling off spatially Gaussian 
with a sigma of 2 Mpc. Finally, we incorporate 
a fourth component for the Fornax cluster, cen- 
tered at (-1.9,-15.0, -13.4) Mpc, with ampli- 
tude 235kms _1 and again a spatial falloff with 
sigma of 2 Mpc. Examples: 

x , y , z v x , Vy , v z v r c v 

-3,18, -2 -590.0,1101.7, -298.6 1209 532 
-35,15,-15 -3551.5, 996.3,-2118.3 4179 272 
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B. Appendix B. The zero point 

We find it difficult to choose a zero point for 
SBF using the various Cepheid distances which 
are available because we are faced with two op- 
tions, each with virtues and faults. SBF is best 
measured in elliptical and SO galaxies, and there 
are no direct Cepheid distance to such galaxies. 
Worse, there are presently no distances to spiral 
galaxies which appear to be tidally interacting 
with E or SO galaxies (such as NGC 4647 and 
NGC 4649). We therefore are faced with the un- 
pleasant choice of associating Cepheid distances 
to SBF galaxies via group associations, or else 
measuring SBF directly in spiral bulges where 
Cepheid distances exist. 

B.l. Option 1. Calibrating via Groups 

For a group comparison we refer to Ferrarese 
et al. (1999) for Cepheid distances, but we reject 
the Coma groups because of the extreme confu- 
sion in that region, the Cen A group because it 
has a very large depth inferred from its angular 
extent across the sky, and the NGC 7331 group 
because it is really a galaxy-galaxy comparison 
with all the inherent worries described below. We 
also, for the sake of argument, do not agree with 
Ferrarese et al. in the inclusion of NGC 1425 
as part of Fornax. It is statistically inconsistent 
with the other two Cepheid galaxies, NGC 1326A 
and NGC 1365; it is halfway between the Fornax 
and the Eridanus clusters, and its Cepheid dis- 
tance is 0.40 mag more distant than the other 
two, very close to the 0.46 mag that we find for 
the mean distance between Fornax and Eridanus. 
We also concur with Ferrarese et al. in rejecting 
NGC 4639 from the Virgo cluster because it is 
0.77 mag more distant. 



Table Bl lists the mean mj we find for these 
groups, adjusted to a fiducial color of (V— I)o = 
1.15, the distance modulus from Ferrarese et al., 
and the inferred SBF zero point, the absolute Mj 
at the fiducial color. 

These groups paint a very consistent story: 



the zero point of the SBF calibration is Mj = 
—1.61 ±0.03. The big advantage to using a group 
calibration is that we have many, very good SBF 
measurements in groups, as is reflected in the 
very small error bars. 

The biggest disadvantage is that the associa- 
tion of spirals with ellipticals is shaky at best. As 
pointed out by Jacoby et al. (1992) the spirals in 
groups seem to lie outside of an elliptical core, as 
one might expect from the fragility of spirals and 
the harsh environment in the center of a group. 
However, in this case we are finding a sheet of five 
Cepheid bearing spirals at precisely the distance 
we assign to the core of the Virgo cluster, and 
two spirals at precisely the same distance as the 
core of Fornax. Had we included NGC 1425 as 
part of Fornax and NGC 4639 as part of Virgo, 
we might have expected the elliptical core to lie 
somewhere in between, leading to a brighter zero 
point by 0.1-0.2 mag. 

B.2. Option 2. Calibrating via Galaxies 

We again refer to Ferrarese et al. for Cepheid 
distances to galaxies where we have measured an 
SBF magnitude in the bulge. Table |B^ lists the 
mean mj we find for these six galaxies, adjusted 
to a fiducial color of (V— J)q = 1.15, the distance 
modulus from Ferrarese et al., and the inferred 
SBF zero point, the absolute Mj at the fiducial 
color. 

These galaxies also are self consistent: the 
zero point of the SBF calibration comes out as 
Mj = —1.80 ± 0.08. Given the wide range in 
the errors in the SBF measurement, it is also in- 
teresting to look at the median: Mi = —1.735. 
Unfortunately the galaxy average is inconsistent 
with the group calibration - they differ by 2.2 
sigma which makes us very reluctant simply to 
take the group-galaxy comparison average. 

The advantages of using the very same galax- 
ies to carry out the calibration is that there is no 
doubt about the spatial relationship of the stars 
contributing to the SBF mean and the Cepheids. 
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Table Bl 
Zero Points from Galaxy Groups 



Group 


{mi) 


(m-M) 


{Mi) 


M31 


22.74 ±0.04 


24.44 ±0.10 


-1.70 ±0.11 


N1023 


28.26 ± 0.08 


29.84 ± 0.08 


-1.58 ±0.11 


Fornax 


29.78 ± 0.03 


31.41 ±0.06 


-1.63 ±0.07 


M81 


26.27 ±0.09 


27.80 ± 0.08 


-1.53 ±0.12 


Leol 


28.48 ± 0.04 


30.08 ± 0.06 


-1.60 ±0.07 


Virgo 


29.41 ± 0.02 


31.03 ±0.06 


-1.62 ±0.06 



Table B2 

Zero Points from Individual Spirals 



Galaxy {mi) {m-M) {Mi) 



N0224 
N3031 
N3368 
N4548 
N4725 
N7331 



22.67 ±0.06 
26.21 ±0.25 
28.34 ±0.21 

29.68 ± 0.54 
28.87 ±0.34 
28.85 ±0.16 



24.44 ±0.10 
27.80 ± 0.08 
30.20 ±0.10 
31.04 ±0.08 
30.57 ±0.08 
30.89 ±0.10 



-1.77 ±0.12 
-1.59 ±0.26 
-1.86 ±0.23 
-1.36 ±0.55 
-1.70 ±0.35 
-2.04 ±0.19 
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The disadvantage is that the SBF measurement 
is very hard, as is evidenced by the large error 
bars. We can only work on the side of the bulge 
which is in front of the plane of the disk; we 
must contend with dust in the bulge; we are af- 
fected by a portion of the disk which underlays 
the bulge we are analyzing; we encounter strong 
color gradients; and we always worry that some- 
how SBF might be different in spiral bulges than 
in ellipticals (despite the evidence to the con- 
trary presented in SBF-I) . We can try to improve 
our SBF measurements, but for some galaxies 
(e.g., NGC 4414) even HST resolution only re- 
veals more and more dust threaded throughout 
the bulge. We do our best to make an accurate, 
unbiased measurement of SBF in spiral bulges, 
but we suspect that if there is a bias, it is prob- 
ably in the sense of assigning too bright an SBF 
magnitude. 

Given the probability that the ellipticals and 
spirals in these groups may be systematically sep- 
arated in space, we feel that we must use the 
galaxy-galaxy calibration, although we are hedg- 
ing our choice slightly by using the galaxy me- 
dian. This is a situation which must be improved 
by further observation. 
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